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' Abstract. 

Q ' We derive general constraints on the rehc abundances of a long-hved particle which 

. mainly decays into a neutrino (and something else) at cosmological time scales. Such an 

I exotic particle may show up in various particle-physics models based on physics beyond 

T-H ■ the standard model. The constraints are obtained from big-bang nucleosynthesis, 

, cosmic microwave background and diffuse neutrino and photon fluxes, depending on 

' the lifetime and the electromagnetic and hadronic branching ratios. 

(N ■ 
in ; 

O I 1. Introduction 

o . 

• In modern cosmology, success of the big-bang nucleosynthesis (BBN) and existence of 
' the cosmic microwave background (CMB) are important facts that support the standard 
' big-bang model. Prediction of the standard BBN scenario is in a reasonable agreement 
' with the observations and the COBE observations P showed the perfect blackbody of 
CMB spectrum. On the other hand, in particle physics, if we consider physics beyond 
the standard model, there exist various new particles some of which have long lifetimes 
and decay during or after BBN. Examples of such long-lived particles include gravitino 
and moduli predicted in the framework of supersymmetry (SUSY) and string theories. 

BBN and CMB are useful probes to exotic particles predicted in physics beyond 
standard model. In fact, the prediction of BBN changes significantly if there exists 
an exotic massive particle with long lifetime. (Hereafter, we call such a particle X.) 
When the lifetime of X is longer than about 1 sec, the decay of X may induce 
electromagnetic and hadronic showers, which lead to photo- and hadro-dissociation 
of ^He and subsequent non-thermal production of other light elements (D, ^He, ^Li, 
and ''Li). Such processes may significantly change the prediction of the standard BBN 
scenario and, consequently, resultant abundances of light elements may conflict with 
observations. Furthermore, the electromagnetic energy injection causes distortion of 
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the CMB spectrum. Since the observation P shows that this distortion is quite small, 
we can constrain the abundance of X. Finally, if the lifetime is very long, the spectrum 
of neutrinos and photons produced by the decay of X are not thermalized and may be 
directly observed. 

The effects of the long-lived particles on BBN were well studied for radiative 
decay [21 El 111 E] and for hadronic decay [6l [71 |8l [9l [TOl HH [12] and stringent constraints 
on the abundance and lifetime of X were obtained. However, when X mainly decays into 
neutrinos, it is expected that the constraints become much weaker because of weakness 
of interactions between neutrinos and other standard model particles. The specific 
case where a sneutrino mainly decays into a gravitino and a neutrino has been already 
discussed in [TH |15l [16]. Such a scenario is realized when the sneutrino is the next 
lightest superparticle (NLSP) while the gravitino is the lightest superparticle (LSP)I3 
In such a case, interaction of X (i.e., sneutrino) is well-known, and it is found that BBN 
provides the most stringent constraint. 

However, there are other possibilities of having long-lived massive particles which 
dominantly decay into neutrinos. For example, if the LSP is axino (a) and the NLSP 
is sneutrino, dominant decay process of the sneutrino (= X) is v v + a. Decay rate 
of this process depends on the properties of axion supermultiplet. Thus, in general, 
properties of X (i.e., lifetime, hadronic branching ratio, and so on) is model-dependent. 
Consequently, the most stringent bound may not be from BBN. For example, high 
energy neutrinos emitted in the X particle decay was considered in |13j where the 
upper bounds on the X abundance were obtained from nucleon-decay detectors and 
Fly's Eye air shower array. 

In this paper, we derive general cosmological constraints on scenarios in which 
there exists a long-lived massive particle which dominantly decays into a neutrino (and 
something else). We treat the lifetime and hadronic branching ratio of the long-lived 
particle as free parameters. In this case, in fact, the constraint from the main decay 
mode is quite weak, and other subdominant decay channel which contain electronic and 
hadronic particles may be important [16]. In our analysis, we take into account both of 
these decay channels and discuss various cosmological constraints. 

Organization of this paper is as foUowis. In Section [2, we discuss BBN constraints. 
In Section [3l constraints from CMB is considered. Then in Section [H we study 
constraints from diffuse neutrino and photon. Section [5] is devoted for conclusion. 

2. BBN Constraints 

Before going into the main subject of this section, which is the constraints from the 
BBN, let us first summarize the properties of the long-lived heavy particle X which 
mainly decays into a neutrino (and some other weakly interacting particle). In our 

X For the case that the gravitino is the the NLSP and that sneutrino is the LSP, see also [14j . 
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study, we assume that X mainly decays as 

X ^ z/ + F, (1) 

where Y is an invisible particle which is very weakly interacting so that it does not 
cause any subsequent scattering with background particles. One of the well-motivated 
examples is the case where the sneutrino is the NLSP while gravitino (or axino) is 
the LSP. (Then, X is the sneutrino and Y is the gravitino or axino.) In addition, for 
concreteness, we assume that the final-state neutrino is electron neutrino. (We note 
here that we have checked that the constraints on the properties of X are not sensitive 
to the flavor of the final-state neutrino.) 

Even though the dominant decay mode is the two-body process, one should keep 
in mind that decay chennels with three- and/or four-body final state should also exist 
since the neutrino as well as X and/or Y couple to Z- and VT-bosons. The emitted 
(real or virtual) weak bosons subsequently decay into quarks and leptons. With this 
type of three- and/or four-body decay processes, energetic quarks and charged leptons 
are produced. 

If the decays of X occur during or after BBN, the standard-model particles emitted 
in the decay can affect the abundances of primordial light elements. First the high energy 
neutrinos emitted in the main decay mode (two-body decay) scatter off the background 
leptons and produce charged leptons (e^, /i^) and charged pions. The former induce 
electromagnetic showers which destroy light elements, and the latter change the n-p 
ratio through nucleon and pion interactions. Second X decays into electromagnetic and 
hadronic particles via the three- and/or four-body decay modes with small branching 
ratio Bx- Such processes directly induce electromagnetic and hadronic showers and 
change the abundances of light elements [10] . 

2.1. Two-body Decay 

First, we discuss effects of the dominant decay process X ^ Y + Ug. We presume that 
Y produced in the decay is a very weakly interacting particle, and that it is irrelevant 
for BBN. Neutrino, however, may affect abundances of light elements. The emitted 
energetic neutrinos scatter off background leptons via weak interaction and several kinds 
of particles may be pair-produced. 

First, charged leptons may be produced via the following processesl§| 



J^i + i^i,BG ^ e~ + e~^ (2) 

+ i^i,BG -> /U" /i+ (3) 

+ i^e,BG ^ yU" + e+ (4) 

+ i^^,,BG ^ e" -h /i+ (5) 



§ Neutrinos may also scatter off the background electron and positron. However, since the photo- 
dissociation processes become important when the cosmic temperature becomes much lower than 
1 MeV. At such temperature, number densities of electron and positron are exremely suppressed, 
and hence the scattering processes with electron and positron are irrelevant for the production of 
charged leptons. 
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where i = e, /z, r is flavor index, and the subscript "BG" is for background particles. 
The muons emitted in the above processes quickly decay into electrons and neutrinos. 
Thus, the above processes produce energetic electrons and positrons which cause 
electromagnetic cascade. Energetic photons in the cascade induce photo-dissociation 
processes of light elements. Effects of these processes have been already studied in |14j . 

Another possible effect is due to the production of pion pairs. High energy neutrinos 
scatter off the background neutrinos and electrons (positrons) and produce pions as 

I'i + i^i,BG ^ TT" + 7r+, (6) 
^. + <BG -^^° + ^^- (7) 

The nucleus-pion interaction rate is ~ lO^sec"^ x (T/MeV)^ which is larger than the 
decay rate of the charged pion (~ 4 x lO^sec"^) for T ~ IMeV. Therefore, the charged 
pions produced at T ~ IMeV scatter off the background nuclei and change protons 
(neutrons) into neutrons (protons) via 

vr" +p ^ n + 7r°, n + 7, (8) 
vr^ + n — > p + 7r°, p + 7. (9) 

Consequently, n/p ratio is increased, resulting in more ^He. Notice that, because of very 
short lifetime, the neutral pions decay before they scatter off the background nuclei and 
hence they are harmless. 

In order to estimate effects of the high-energy-neutrino induced processes, we have 
numerically solved the Boltzmann equation describing the time evolution of the high 
energy neutrino spectrum taking into account of all the processes above. Details of our 



calculation are given in Appendix A 



2.2. Three and Four-body Decay 

Even though the branching ratio for three and four-body decay processes are much 
smaller than 1, such decay processes are very important since colored and charged 
particles are directly emitted from these decay processes. Energetic colored and 
charged particles may significantly change the prediction of the standard BEN scenario. 
Effects of these particles are classified into three categories: photo-dissociations, hadro- 
dissociations and p n conversion. 

In order to study the effects of three- and four-body decay processes, it is important 
to obtain the spectra of quarks and leptons emitted by the decay of X, which depends on 
the model. In our analysis, we use, up to normalization, those obtained in the case where 
X is the sneutrino and Y is the gravitino; in such a case, the sneutrino may decay into 
the gravitino, neutrino (or charged lepton), and Z^*'^ (or W^*^), and the produced (real 
or virtual) weak bosons subsequently decay into quark- or lepton-pair. (Here, Z* and 
W* denote virtual weak bosons.) In order to perform our analysis as model- independent 
as possible, we treat that the branching ratio for three- and four-body processes is a 
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Figure 1. Spectrum of e+ + e (solid line) and photon (dashed line) with mx 
lO^GeV. 



free parameter: we define branching ratio Bx as 

r(X 3 body) + r(X ^ 4 body) 



B 



X 



r(x all) • (1°' 

where F is decay width. For the study of the effects of photo-dissociation processes and 
CMB spectral distortion, we also calculate averaged "visible energy" emitted from one 
X: 

-Evis = Bx{E^is), (11) 

where (-Evis) are averaged energy carried away by charged particles and photons in 
three and four-decay modes. The energy distributions of photons, neutrinos, leptons 
and nucleons produced by the decay of X are calculated by means of Monte Carlo 
simulations with the PYTHIA package [17]. In Figs. [Hill we plot the spectra of photons 
and leptons (e+ + e~), respectively. From these distributions we calculate the averaged 
energy and obtain (E^is) = 25.3, 146, 821 and 5630 GeV for mx = 10^, 10^, 10^ and 
10^ GeV, respectively. Here, mx is the mass of X. E^^g is much smaller than mx since 
we are interested in the case where Bx <C 1. 

Photo-dissociation processes are induced by energetic photons in electromagnetic 
shower which is caused by charged particles and/or photons emitted from X. With 
given background temperature, the distribution function of energetic photons depends 
on total amount of energy injected by particles with electromagnetic interaction, and 
is insensitive to the shape of energy spectrum of primary particles. Thus, once i^vis 




Figure 3. Same as FigH] except mx = lO^GeV. 
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Figure 4. Same as Fig[T] except mx — lO^GeV. 



is obtained, the energy distribution of energetic photons in the electromagnetic shower 
can be obtained. Then, photo-dissociation rates are obtained by convoluting energy 
distribution function and cross sections of photo-dissociation reactions. For details of 
our treatment of photo-dissociation processes, see |3j. 

For the study of hadro-dissociation processes, it is necessary to obtain energy 
distributions of (primary) hadrons which are produced after the hadronization of quarks 
emitted from X. We have calculated the spectra of p and n by using PYTHIA. 
These hadrons cause hadronic shower and induced hadro-dissociation processes. In 
our analysis, in addition, we have also calculated the number of charged pions produced 
by the decay of X. Such charged pions, protons and neutrons become the source of 
p ^ n conversion process, which changes the number of ^He [71 [9] {JJ] Once the spectra 
of hadrons are obtained, effects of hadro-dissociation and p ^ n conversion are studied 
with the procedure given in [lOj . 




2.3. Numerical Results 

In our analysis, we have followed the evolutions of the number densities of the light 
elements. For this purpose, we have modified the Kawano code [H] including photo- and 
hadro-dissociation processes. As observational constraints on the primordial abundances 
of light elements, we adopt those used in [16] except for Yp and (n3He/'"'D)p: 



II Here we have neglected the effects of Kaons to the p ^ n conversion, according to the discussion 
in [10]. 





(2.82 ±0.26) X 10~^ 



(12) 
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(n3He/nD)p < 0.83 + 0.27, (13) 

Fp = 0.2516 ± 0.0040, (14) 

logio(n7Li/nH)p = -9.63 ±0.06 ±0.3. (15) 

(n6Li/n7Li)p < 0.046 ±0.022 ±0.084, (16) 

where the subscript "p" is for primordial value (just after BBN), and Fp is the primordial 
mass fraction of ^He. For the center value of Yp, we have adopted the value reported 
in [19] in which the authors used new data of Hel emissivitiesj^l and conservatively 
added a larger error (= 0.0040) as discussed in [21]. For (n3He/''^D)p5 "we have adopted 
most newly-reported values of D and ^He abundances observed in protosolar clouds [22], 
(n3He/nH)psc = (1.66 ±0.06) x 10"^ and (nD/nH)psc = (2.00 ±0.35) x 10-^ where the 
subscript "PSC" means a value in the protosolar cloud. (For the importance of the 
upper bounds on (n3He/''^D)p, see [ZHlHOj.) 

We parameterize the primordial abundance of X by yield variable Yx which is 
defined as the ratio of number density and total entropy density at (t <^ rx); 

Yx - PI , (17) 

where tx is the lifetime of X. If Yx is too large, abundances of light elements are too 
much affected to be consistent with the observations. Thus, we can derive upper bound 
on Yx- 

In Figs. [SlfTOl we show BBN constraints on Tx and nixYx for Bx = 10^'^ and 
Bx = 10~^; we found that constraints with Bx = is almost same as that with 
Bx = 10-6. 

As one can see from these figures, the most important constraint comes from 
overproduction of ^He when tx ^ lO^sec. Since protons are more abundant than 
neutrons, a significant amount of proton may be converted to neutron through 
nucleus-pion interaction (and nucleon-nuclen interactions for relatively large Bx) and 
consequently "^He is overproduced. When lO^sec ^ tx < lO^sec, the background "^He 
(which we call obg) is effectively dissociated by the energetic hadrons produced in the 
hadronic shower. In this case, overproduction of D may occur as a result of hadro- 
dissociation of obg- In addition, energetic T and ^He are also produced and they 
synthesize ^Li through the ^Li via T ± obg ^Li ± n and ^He ± aec — ^Li ± p. 
When Tx ^ 10'' sec, the energetic hadrons are stopped by the scattering processes with 
background electrons, and hence the effects of hadro-dissociation become less efficient 
than those of the photo-dissociation. In particular, the energetic photons produced in 
the electromagnetic shower destroy obg- In this case, overproduction of D and ^He 
occurs as a result of photo-dissociation of aBc- 

The constraints on mxYx depend on mx in a non-trivial way. The contraint 
coming from the photo-dissociation caused by the two-body decay becomes stringent 
as mx increases because neutrinos emitted in two-body decay have higher energy and 

^ See also the other recent value of yp reported in [20] where the authors adopted larger errors (0.0028) 
than that of [HI. 
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Figure 5. BBN constraints on tx vs. mx^x plane. Here, we take mx — lOOGeV 
and Bx = 10^^. 

scatter off the background neutrinos with larger rate. On the other hand, the rates of 
the hadro- dissociation and photo-disocciation caused by the three- and four-body decay 
depend only on E^i^x and hence the constraints on mxYx become slightly milder as 
mx increses. The constraint from n ^ p conversion becomes weaker as mx increases. 
This is because the charged pion production is roughly determined by Yx for the three- 
and four-decay. 

So far, we have shown the results for the case where the primary neutrino is electron- 
type. However, we have checked that the BBN constraints are almost unchanged even 
if X decays into muon- or tau- neutrino (and Y). 

3. CMB Constraints 

CMB also imposes constraints on the decays of X. COBE observations show that 
CMB spectrum is almost perfect blackbody [Ij. Therefore any exotic energy injections 
that cause distortions in the spectrum of CMB are stringently constrained [211 [25l [26] . 
When the photons are emitted before redshift z ~ 10^, they are thermalized by Compton 
scattering, double Compton scattering and bremsstrahlung, and no spectral distortion 
takes place. However, at z < 10^ only Compton scattering is efficient. Since the 
Compton scattering does not change the total number of photons, the resultant spectrum 
becomes a Bose-Einstein distribution with a finite chemical potential /x, regardless of 
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Figure 6. BBN constraints on tx vs. mxYx plane. Here, we take mx = ITeV and 
Bx - 10-3. 
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Figure 7. BBN constraints on tx vs. mxYx plane. Here, we take mx = lOTcV and 
Bx = 10-3. 
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Figure 9. BBN constraints on tx vs. mxYx plane. Here, we take mx — ITeV and 
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Cosmological Constraints on Neutrino Injection 12 




10-1 1 10 lo^io^iono^ionono^io^ioi' 

Xx (sec) 



Figure 10. BBN constraints on tx vs. mxYx plane. Here, we take mx = lOTeV 
and Bx = 10"^. 



the detail of the injection. For the case of massive particle decay, the number density 
of injected photons is negligible compared with that in the background. Therefore, the 
spectral distortions are determined by the fraction of the energy release, Ap^/ p^. Then, 
for small Ap^/ the chemical potential of the photon spectrum is given by 

0.714 p^ 

As we have discussed, there are two types of processes which contribute to the 
electromagnetic energy injection (Ap^) in the present scenario. One is the two- 
body decay process: the high energy primary neutrinos emitted by X scatter off the 
background neutrinos and create charged leptons whose energy is finally converted to 
the energy of radiation. The other is the three- and four-body decay processes by which 
quarks and charged leptons are produced. In our analysis, we have taken into account 
both of these contributions. 



3.1. Three and four-body Decay 

There are two quantities which are required to calculate the electromagnetic energy 
injection: branching ratio Bx and the averaged energy (-Evis)- 

We follow the treatment of [26] for the case of three and four-body decays. Using 
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the fact that X decays exponentially in time with lifetime r^, we obtain 

-Bx, (19) 



Ap^ _ (E^vis) nx 



Pj 2.701T(tofr) 

where T{t) is the CMB temperature and nx is the number density of X before decay. 
Here tefr = [r(l + /3)]^''^ tx for time-temperature relation T oc t~^, where T is the 
gamma function. 

From Eqs. (fTSll and ( fTOl) . we find that the chemical potential is given by 

..4.oox:o'^(^)-(iM)..^. (.0, 

We have assumed here that we are in the radiation dominated epoch where T oc 
Note that, however, photon number changing processes (double Compton scattering 
and bremsstrahlung) become increasingly efficient as the photon frequency decreases. 
This means the spectrum becomes blackbody at low frequencies. The photons with 
low frequencies produced by the photon-number changing processes are transferred to 
higher frequencies by inverse Compton scattering and the chemical potential decreases 
in time. For a low Qbh"^ universe suggested by the BBN [29] and WMAP [30J (where, in 
this paper, Qj, denotes the density parameter of baryon, and h is the Hubble constant 
in units of 100 km/sec/Mpc), double Compton scattering dominates the thermalization 
process. The chemical potential produced a.tt = th is blurred out at an exponential rate 
and the present value is given by [28] 

/io = /i(to) = fiih) exp{-{tBc/thf) (21) 



with 



/ O 1.2 \ 4/5 / -1^ \ 4/5 

toe = 6.81 X lOM 1-- sec. (22) 



Combining Eq. (l20l) with Eq. (1211) . we find that the chemical potential today is given 
by [26] 

/.o ^ 4.00 X 10^ (^Y'eM-itBc/rxf') (^) (23) 
Msec/ V ICeV / 



ViGevy 

For late energy injection {z < 10^), Compton scattering can no longer establish the 
Bose-Einstein distribution. In this case, the spectrum can be described by the Compton 
y-parameter which is defined by 
r T - T 

y = / dt^ n.ar, (24) 

J 

where me and Te are the number density and temperature of electrons and ct is Thomson 
cross section. Then the energy injection is related to y as Ap^/p^ = Ay. Here, we 
define zk as the redshift at which the time scale for energy exchange through Compton 
scatterings is equal to the Hubble time; 
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The spectrum can be described by chemical potential /i for energy injection at ^ > 2;^^ 
and by Compton ^/-parameter for energy injection at 2; < zk- Here we take VLi^h"^ = 
0.0223 [30], < 9 X 10"^ and \y\ < 1.2 x 10"^ as observational limits by COBE f3n[32]. 



3.2. Two-body Decay 

In addition, the effect of charged-lepton productions through scattering of high energy 
neutrinos off background neutrinos (see Eqs. ([2]) - ([5])) should be taken into account. 
The amount of energy which is converted to the background photons is estimated by 

TT = / 77^;7TT77; — r{t,m,^,Tx), (26) 



P7 Jo P7 dt 2 Jq 2.701Trx 

where dEi/dt is the energy density which is converted to charged lepton and finally 
photon through neutrino scattering per unit time and 

\ 2tx J dt ^ ^ 

Notice that r represents the ratio of the radiative energy injection per X decay to the X 
mass. The details of calculation are found in [m [33]. In Figs. [11] and [121 we show the 
time evolution of r. Fig. [TT] shows that the ratio r increases with mx. This is because 
higher energy neutrino has larger cross section for scattering off the background neutrino 
and is also easy to exceed the threshold energy for lepton pair creations. The similar 
logic applies to Fig. [121 In this case, background neutrino energy increases with the 
decrease of tx, so that the ratio r increases. 



3.3. Constraint from CMB 

As mentioned earlier, the spectral distortions are determined by Ap^/ which is the 
sum of Eqs. ( [T9l) and ( l26l) . In Figs. [131 and El we show the upper bounds of mxYx 
taking account of neutrino-neutrino scattering. 

When the lifetime is short and the energy of background neutrino is sufficiently high, 
H is determined by neutrino-neutrino scattering, hence r. Therefore, the constraints 
become severer as mx becomes larger. When lifetime is long, is determined by three 
and four-body decay of X, hence {E^i^) jmx- As already mentioned, {E^i^) jmx becomes 
smaller with larger mx- Then, the constraints become severer as mx becomes smaller. 

So far, we have focused only on photon energy injection. However, the emitted 
ultra-relativistic particles (neutrino and Y) contribute to the total relativistic energy 
and could lead to a more stringent constraint than that from spectral distortion, when 
branching ratio is sufficiently small. Before recombination, the CMB angular power 
spectrum is sensitive to the change of the total relativistic energy through the early 
integrated Sachs- Wolfe effect [M]- In addition, too much relativistic energy affects the 
growth of large scale structure (LSS) since the epoch for the matter-radiation equality 
becomes later. The increase of the total relativistic energy is conventionally described 
by an effective number of light neutrino species /S.Ny. The combined analysis of CMB 
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Figure 11. The time evolution of r with mx — lO^GeV (tliin sohd hne), lO'^GeV 
(thin dashed line), lO'^GeV (thick solid line) and lO^GeV (thick dashed line). We take 
Tx = lO^sec. 




Figure 12. The time evolution of r with tx = lO^sec (solid line), lO^sec (dotted line), 
10"sec (dashed line). We take mx = lO^GeV. 




Figure 13. CMB constraints on mxYx with Bx = 10 ^. From upper to lower 
the hnes represent the upper bound of mxYx when mx = lO^GeV (thin soHd hne), 
lO^GeV (thin dotted hne), lO^GeV (thick sohd hne) and lO^GeV (thick dotted hne). 




Figure 14. Same as Fig. [T3l except Bx = 10 ^. 



and LSS data sets the upper bound on AA^j, as AA^j, < 5.0 [35]. In Fig. [T5l we show the 
constraints from the CMB spectral distortion with mx = 10^ GeV when Bx = 10^^ and 
10~^. In addition, we also show the constraint from total relativistic energy injection. 
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which is independent of Bx- When branching ratio is sufficiently small, Fig. [13] shows 
that the constraint from the total relativistic energy injection provides severer constraint 
in a wide range of lifetime. 
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Figure 15. Solid and dotted lines represent the CMB constraints on mxYx with 
rax = 1 TeV when Bx = 10^^ and 10~^ respectively. Dashed line represents the 
constraints from increases of the total relativistic energy. 



4. Diffuse neutrino and photon Constraints 
4-1- Diffuse Neutrino Flux 

When neutrino injection takes place very late, the emitted neutrinos may produce 
an observable peak in the diffuse neutrino spectrum. The present differential flux of 
neutrinos is given by [Ml ET] 

d^. If, 1 YxSq dN, 

-TT^ = — / dz———— exp — t/Tx -TT^, 28 

where Sq is the present entropy density and h{z) = [(1 + z)^Qm + ^^a]^^^ with fi^ 
and f^A being the density parameters of non-relativistic matter and dark energy, 
respectively. Here neutrinos produced with energy E are redshifted to the observed 
energy Eq = E/{1 + z) and E is just half of the mass of X. The source spectrum 
dN,y/ dEo is given by 

^^ = /.5(E-Eo(l + .)), (29) 

where fi is the fraction of the neutrino spades i = e, ij,,t emitted by the decay of X. 
(Notice that = 1.) We have assumed that only electron neutrinos are produced 
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in the decay process. However, we should take the effects of neutrino oscillations into 
account. Neutrino oscillations can be described by six parameters: two independent 
mass differences (AmfgjAmlg), three mixing angles (6'i2, 6*23, 6*13) and a CP-violating 
phase 5. A mixing angle ^23 is ~ 45° from atmospheric neutrino experiments [38]. A 
mixing angle 9i2 is determined by solar neutrino experiments as 612 — 34° |39]. CHOOZ 
experiment presented a mixing angle 61^ < 12° [40]. In our case, neutrino traveling 
distance is very long and mass differences are irrelevant. CP- violating phase 6 enters 
the mixing matrix only in combination with sin6'i3. In a reasonable approximation, 
/e ~ 0.6 and 0.2. 

The present atmospheric neutrino u^^ + z/^ data gives the upper bound of the 
differential flux of neutrinos. The atmospheric neutrino has been observed by 

Super-Kamiokande |38j and AMANDA ^1]. For observational flux of energy range 
0.3 — 1.0 X 10'^ GeV we adopt the result in [12] where the atmospheric neutrino fluxes 
are estimated from the data on atmospheric neutrino event rates measured by the Super- 
Kamiokande experiment. For higher ernergy range 1.3 x 10^ — 3.0 x 10^ GeV we use 
the atmospheric neutrino spectrum derived from AMANDA. In Fig. [16], we show the 
atmospheric neutrino fluxes from the Super-Kamiokande and AMANDA experiments 
as well as the diffuse neutrino fluxes from the X decay. 

Futhermore, the diffuse neutrino flux is also constrained from null detection of the 
relic supernova z/g flux by Super-Kamiokande. In [13], the upperbound on z/g flux is 
obtained as < 1.2cm~^s~^ above threshold of Ey > 19.3MeV. 

We require that the neutrino flux from X decay should not exceed the observed 
atmospheric + flux and upperlimit of the relic supernona z/g flux, which leads to the 
constraints on the abundance of X as shown in Fig. [171 In the figure we also show the 
constraint from relic supernova search only. When lifetime is short, the constraint is very 
weak since the neutrinos get redshifted until the present and their energy becomes lower 
than the 20 MeV. For intermediate lifetime, the constraints on the abundances of X are 
determined by the relic supernova search and they are in proportion to nix because 
the diffuse z/g flux is determined by the number of injected neutrinos above threshold 
energy, hence Yx- When lifetime is long, however, the constraints are determined by 
atmospheric neutrino fluxes. Since the neutrino flux has a peak at ~ mx/(l + za) {zd. 
redshift at t = tx), the maximum differential neutrino flux is propotional to Yx/mx- 
On the other hand, the observed differential neutrino flux is roughly propotional to E^^ 
and hence at ~ mx- Therefore, the constraint on mxYx depends on mx as 
mxTn'^ /m'^ ~ m^. This means the constraints become severe with larger mass as 
shown in Fig. [171 Constraints for lifetimes longer than present time to scale by a factor 
to/Tx relative to the constraints at to = tx- 

4-2. Diffuse Photon Flux 

High energy photons and electrons (positrons) produced in the three and four-body 
decay may be observed as diffuse gamma rays when the decay takes place after the 
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Figure 16. Atmospheric neutrino flux. Thin solid lines represent Icr range of the 
atmospheric neutrino fluxes [42]. The point data are from AMANDA. Thick solid and 
dotted lines represent diffuse neutrino signal with mx — lO'^GeV and Yx — 2.53 x 10"^'' 
(thick sohd hnc) and lO^GeV and Yx = 3.37 x IQ-"^^ {tUck dotted line). The lifetime 
is Tx — lO^^sec. 




Figure 17. Constraints from diffuse neutrino flux. From upper to lower the lines 
represent the upper bound of X abundance when mx — lO^GeV (solid line), lO'^GeV 
(dashed line), lO^'GeV (dotted line) and lO^GeV (dot-dashed line). Thin lines are 
constraints only from Pg. 
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recombination epoch. In calculating the diffuse photon spectrum, we must consider 
primary photon spectrum which is not monoenergetic. In addition, we should also take 
account of the absorption of gamma rays along the line of sight. 

Since photons are produced through three- and four-body decays, their spectrum 
is not monochromatic unlike neutrino. The energy distributions of photons, neutrinos, 
leptons and nucleons produced by the three- and four-body decay of X are shown 
in Figs. [Hill The energy of produced electrons and positrons is transferred to the 
background photons through inverse Compton process. In calculating the photon flux, 
we have taken into account the photons produced by the inverse Compton process as 
well as those from the cascade decay chain induced by the three- and four-body decay of 
X. The latter effect becomes more important for high energy photons. However, since 
the inverse Compton process produces many soft photons, the formar process becomes 
more significant for low energy photons. The details of calculation of inverse Compton 
process are given in Appendix B Nucleons are also produce in the decay and they 



produce photons through inverse Compton process or 13 decay. However, we neglect this 
effect since the number density of produced nucleons is sufficiently small. In Fig. [181 we 
plot photon flux for mx = iC^GeV. 
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Figure 18. Photon spectrum versus photon energy for mx = lO^'GeV. Sohd hne 
represents photon produced through inverse Compton process at 1 + z = 1 and dotted 
Hne represents at 1 + z = 100. Dashed hne represents photon produced by direct decay 
of X. 



High energy photons injected in the universe, in general, scatter through various 
processes; photon pair production, photon-photon scattering and pair production in 
matter. At the eraly epochs absorption and scattering due to the background photons 
are important, whereas significant absorption by diffuse IR-UV photons emitted from 
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galaxies takes place at later epochs. 

First let us consider the radiative processes due to the background photons. In 
the present case, the relevant processes are photon-photon scattering and photon pair 
creation [33]. Then, photons degrade their energy by producing electron-positron pairs 
or dividing their energy with the background photons. The photon spectrum was 
calculated in detail in |45j. According to [15] the present differential flux of photons 
is given by 



where Ci{E, z) is the number of photons per unit energy produced both directly and by 
inverse Compton at z for one X decay, and 700) is the redshift at which the optical 
depth of high enrergy photons becomes 1. In addition, £s{Ei, E2, z) is the number of 
photons per unit energy in the spectrum when photons with energy Ei are produced by 
scattering of photons with energy E2 at redshift z. If there is no scattering, Cs{,Ei, E2, z) 
becomes 5{Ei — E2). (The concrete expression of Cs{Ei, E2, z) is found in [35].) 

In addition, 7 rays with GeV to TeV energies are absorbed via electron-positron 
pair production on diffuse background IR-UV photons which have been emitted by 
galaxies [lEl HTl HE]. In this paper, we adopt the result of [3H] which calculated the 
optical depth tir of the universe for 7 rays having energies from 4 GeV to 100 TeV 
at redshifts from to 5. (See Fig. 8 in [3H]-) Then the resultant photon spectrum is 
given by Eq. (l30l) multiplied by e""^^. We neglect the secondary soft photons produced 
via electron-positron pair production on diffuse background IR-UV photons, and only 
consider the attenuation of high energy photons. 

In Fig. UHl we show the diffuse photon flux from the COMPTEL [39] and EGRET 
[5U] observations and decay of X, from which we obtain the upper limit on the abundance 
of X as shown in Fig.[20l Compared with the limit from the neutrino flux, the constraints 
from photon flux are almost insensitive to mx- This can be understood as follows. The 
differential photon flux at the peak energy is roughly proportional to Yx / ttlx from the 
same reasoning as the neutrino flux in Sec. 14.11 while the observed one is proportional 
to E'"^ oc m]^^. Thus, the constraint on Yx depends on mx as rn^mx ~ ""^x^j which 
means that the limit on mxYx is almost insensitive to mx- Fig. [20] shows that the 
constraints become less stringent with larger mx when lifetime is long. There are two 
reasons for this. One is EGRET had observed up to lOOGeV. When mx is large, the 
present photon energy at which the flux becomes maximum may exceed the energy 
range of EGRET observation. The other is higher energy photons are more effectively 
absorbed by diffuse background photons. 




(30) 
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Figure 19. Diffuse plioton flux. Tlic point data are from COMPTEL and EGRET. 
Solid and dotted lines represent diffuse photon signal for mx — 10^ GeV and 
BxYx = 8.86 X 10-2°(solid line) and 10^ GeV and BxYx = 2.00 x 10-2i(dotted 
line). The lifetime is tx — 10^^ sec. 




Figure 20. Constraints from diffuse photon flux. From upper to lower the lines 
represent the upper bound of X abundance when rax = 10^ GeV (thin solid line), 
10^ GeV (thin dotted hne), 10"* GeV (thick solid hue) and 10^ GeV (thick dotted 
line). 
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5. Conclusion 

In this paper, we have considered the long-lived massive particle X which mainly decays 
into a neutrino and an invisible particle, and have investigated the cosmological and 
astrophysical constraints on the high-energy neutrino and photon injection due to decay 
of X-particle. We have shown that the BBN, CMB, diffuse neutrino fluxes and diffuse 
gamma rays provide stringent constraints on the abundance of the decaying particle X. 
We summarize the constraints in Figs. EH and [221 




Figure 21. The constraints on the relic abundance of X from various observations 
with mx = lOOGeV and Bx = 10"^. 
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Appendix A. Bottzmann equation 

In order to investigate effects of photo-dissociation processes and p ^ n conversion 
processes, we have to calculate a photon spectrum. In this paper, a source of high 
energy photon is charged leptons and pions which are produced through scattering 
of high energy neutrinos off background leptons. Therefore, we determine the time 
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Figure 22. The constraints on the rehc abundance of X from various observations 
with mx = lOTeV and Bx = 10"^. 



evolution of the distribution function of high energy neutrinos in order to investigate 
the photon spectrum. In this appendix, we write down Boltzmann equations which 
determines the high energy neutrino spectrum. Our notation is the same as |33j . 

The high energy neutrinos (z/) produced in X decay scatter off the thermal neutrino 
(z/fe) in the background by the following processes; 

(A.l) 



Vi + z/, 
Vi + z/, 

Z/j + V. 
Vi + Z/j 

Vi + Vi 
Vi + Vi 



3 



Vi + Vi 
Vi + Vi 



Vj + V. 



01 



Vi + z/. 



J' 



Vi + z/. 
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(A.2) 
(A.3) 
(A.4) 
(A.5) 

e" + e+, (A.6) 
+ (A.7) 

where index i and j represent e, /i and r with i ^ j. All the amplitude squared 
in these reactions take the following form: 

\M\^ = 32Gl[a{pp'f + bipqf + c{pq'f + dm\pp% (A.8) 

where Gp — 1-17 x 10"^GeV"^ is the Fermi constant, the coefficients a - d depend on 
the individual reaction, p and p' are the initial momenta of high energy neutrino and 
background neutrino, q and q' are the final momenta, and m represents the mass of the 
fermion the in final state. Coefficients for each processes are given in Table lAll 

First, let us consider the neutrino scattering processes Eqs. (1A.1I) - (1A.5I) . Here, we 
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Process 



Vi + Vi 
Vi + 

Vj + z/j- 

Z/i + Vj 
Vi + 

+ c 

/7 + 



2 


1 















9 
1 

1 



{Cv-Ca? (Cy + CA + 2)2 (Cv 



[Cy — 



[Cy + Ca? 









CA)(Cy + 



C^A + 2) 



Table Al. Coefficients a - d for each processes. Index i and j (with « ^ represent 
the generation, Vi \y is the background neutrino of i-th generation, if is the charged 
lepton of i-th generation (in our case, or Cy and Ga are defined as follows: 

Cy = —0.5 + 2sin0^ and Ca = —0.5. Here, 9^ is the Weinberg angle. 



define Ei, as the energy of initial liigli energy neutrino and E'^ as tlie energy of neutrino 
in final state. We also write the distribution function of the background neutrino as 

^2 1 



ME.) 



(A.9) 



27r2exp(E,/T,) + 1' 

where E^, is the energy of background neutrino and T^, is the neutrino temperature. We 
describe the contribution to the time derivative of the neutrino distribution function. 
When Ey.El ^ E^,, the increases of distribution function due to scattering is written 
as 



dt 



4 

3^ 



POO 1 

Gl / dE, — [aEl + b{E, - Elf + cE'l]U{E, 
Je' K 



dEyE^f^{Ey 



(A.IO) 



On the other hand, decreases of neutrino distribution function is written as 



dUE, 



dt 



--f{E 



' - I - - 
dEy^fi,{E, 
Et 



dssa{s) 

3^.3-,, dEMuiE,), (A.ii; 

where a{s) is the total cross section obtained from the amplitude Eq. (lA.Sp . Notice that 
the condition for the neutrino number conservation is realized; 



4 ^2 
on 



-b + -c 



dt 



dt 



(A.12) 



unless the effects of inelastic channels (z/ + z/ e+ + e~,/i^ + /i~) are taken into account. 

Effects of thecharged lepton pair creation process can be taken into account in 
the same way, and the contribution to the time derivative of the neutrino distribution 
function is given by 



dt 



11 /■°° - 1 - - 



dssa(s) 



Cosmological Constraints on Neutrino Injection 



26 



El- 



(A.13) 



with 



4 4m^ 



EyEi, 



1/2 



- 4m^ In 
/i = 2 f 4 - 



E,E,(SElEl - 2m^E,E, - 3m^) 
2{4 - (4mV^.^.)}'/'^.^. + 4E,E, - 2m2 



2m^ 



(A.14) 



4m^ 

E^Ey 



1/2 



E,E,{2E,E, - 



— 2m In 



2{4 - (4mVE,E,)}^/^E,E, + 4E,^, - 2m^ 
2m2 



(A.15) 



Coefficients for the charged lepton production processes are given in Table lAll 
In addition, there are processes for neutrino scattering as follows: 



+ ^e,b + e+. 

The amplitudes squared in these reactions take the form given by 

1 



I A/fP 



128G'^(pg') 
128G|(pg') 



(A.16) 
(A.17) 

(A.18) 
(A.19) 



Effects of these process can be taken into account in the same way. However, it is 
somewhat complicated to calculate Boltzmann equations for these processes due to the 
mass difference between muon and electron. The contributions to the time derivative of 
the neutrino distribution function are given by 



dt 



fuAEu 



dt 
24^ 



oo 



poo 

EJ^E,) / dE,ME,)h. (A.20) 
Jo 



with 



/3 = - [-5x^ - 2xy - 5y^ - 20{x + y) + 64] 



+ {x^ - 3x'^y - 3xy'^ + y^) In 



4 — X — y + 2z 
2Jxy 



-\- \x — y\^ In 



4(x + y) — (x — yY — 2z\x — y\ 



xy 



(A.21) 



where x = m'j^/E^E^, y = m?^/ E^E^ and z = \A — 2{x + y) + \{x — y)^] ^. 



Cosmological Constraints on Neutrino Injection 



27 



Process 


a 


b 


c 




d 


I'e + ^ Ve + e~ 


{Cv + Ca + 2)2 





{Cv - Ca? 


-{Cv 


-Ca){Cv + Ca + 2) 


z/e + Cfl" ^ z/g + e+ 


(C^y - Ca? 





{Cv + Ca + 2f 


-{Cv 


-Ca){Cv + Ca + 2) 


Vi + el Vi + 


{Cv + C^)' 





{Cv - Ca? 




-{cl - CD 


Vi + e'^ ^ Vi + e+ 


{Cv - CaT 





{Cv + Ca? 




-{c'v - CD 



Table A2. Coefficients a - d for eacli processes. Index i represents /i and r, eb is tfie 
background electron and positron. 



Finally, we consider the effect of charged pion pair creation process. The cross 
section for charged pion pair production is given by 



with 



1 / Am"^ 

(t(z/z7 TT+TT-) = Gl{l-2sm^ew?s[l ^1 \F{s)\\ (A.22) 

127r \ s 



4 

m 

[s — miY + miVi 



where rup and Tp is the mass and decay width of the p meson, respectively [51]. The 
contribution to the time derivative of the neutrino distribution function is given by 
Eq. ( 1A.13I) and the cross section above. 

Next, we turn now to neutrino-electron scattering processes. The high energy 
neutrinos also scatter off the thermal electron (e^) and positron (e^) by the following 
processes; 

z/e + e^ ^ z/e + e", (A.24) 
z/e + e+ ^ z/e + e+, (A.25) 
+ % + e", (A. 26) 

z/i + e+ ^ z/i + e+. (A.27) 

where index i represents /i and r. All the amplitude squared in these reactions take the 
same form as Eq. flA.8|) . In this case, however, p and q are the initial and final momenta 
of neutrino and p' and q' are that of background electron (positron). Coefficients for each 
processes are given in Table IA2[ We write the distribution function of the background 
electron (positron) as 

- - f 1 

fJE^ = At , , . (A.28) 

27r2 exp(Ee/r^) + 1 ^ ' 

where and Pe is the energy and momentum of background electron (positron), 
respectively. Photon temperature is different from neutrino temperature due to 
neutrino decoupling and subsequent electron-positron pair annihilation. From entropy 
conservation, the relation between them is given by 



4 



1/3 



'7 



45 1 p / - p\ 



1/3 



(A.29) 
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The effects of neutrino-electron scattering are only important at early time because of the 
Boltzmann suppression of the distribution function of background electron (positron). 
When El,, El ^ Ee, the increase and decrease of distribution functions are given by 

dMK) 



dt 



-Gl / dE, — [aEl + b{E, - Elf + cE'l]UE, 
Stt Je, E^ 



dt 



dpeEefe{Ee 



1 - 

dpe 



ml 
1 

PeEe 



(A.30) 



fe{Ee) / ds{s - ml)a{s) 



Stt 



Gl(a + -h + -c]EJ,{E, 



dpeEefe{Ee 



1 



ml \ 

aeI) 



(A.3i; 



Notice that the condition for the neutrino number conservation is also realized just as 
neutrino-neutrino scattering. 

We also include the following process: 

+ e~ ^ Ue + fi~ ■ (A. 32) 

The amplitude squared in this reaction is give by 

1 



\M\ 



128GUpp') 



(A.33) 



For simplicity, we neglect the mass difference between muon and electron in this reaction. 
On this assumption, this process is the same form as z/j + e — z/j + e. 

In addition, we consider the effect of pion pair creation process. The cross section 
for pion pair production is given by 



127r 



G'^pS 



Amz 



\F{s) 



(A.34) 



where |F(s)p is defined in Eq. (]A.23|) [51j. The contribution to the time derivative of 
the neutrino distribution function is given by Eq. (]A.3ip with the cross section above. 

Then one can obtain the Boltzmann equations describing the evolution of the 
spectra for the high energy neutrinos; 



dt 



37r 



dEu^ \^in,ijEl + bin,ij{E^ - Elf + Ci^^ijE'l f^^ {E^ 



X 



dE^Eyfy{Ey 







4(^2 

Stt 



out 



^ 2 "out ^ ^'^out 



dEyEyfy{Ey 



4G| 
37r 

oo 



OO -j^ 

dE. 



E?. 



(^nA^u + KnAEv " El) + C-^.-E'^ fuA^u) 



X 







1 77)2 

dPeEME,) [l--^-f. 
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-J^KUM) [al^t, + -/out, + gCL^.j dpeEj^iE,) i^l - 



dt / ,- V dt / ,- -,+ 

+ :^nxS{El-mx/2)Si,e 
Ztx 

+ KH^^0^-2HUK), (A.35) 

where H is the expansion rate of the universe and Sij is a Kronecker delta. The 
coefficients for neutrino-neutrino scattering are given by 

C* = 4, 6^„, = 0, C, = 13, (A.36) 

^in,ii — 6, ^'m,M — 9) Cj„ jj = 11, (A. 37) 

<n,ii = l, ^'rn,^,— !, C,ii = 2, (i^j). (A.38) 

For example, let us derive a^^^. Factor 2 comes from z/jZ/j — >■ i/jZ/j scattering and another 
factor 2 comes from UiUj UiUj scattering. Consequently, the coefficient a^^^ amount to 
4. Other coefficients can be derived in the same manner. The coefficients for neutrino- 
electron scattering are given by 

<ut,e = (Cv + Ca + 2f + {Cv - Ca)\ = 0, = a^„,,„ (A.39) 

<ut,j = {Cv + CAf + {Cv-Ca)\ &Ltj = 0, Ct,j=<„t,„ (A.40) 

"'out,f> "in,i "out,ii ^in,i out,ii \ ) 

where index i represents e, fi and r and j represents fi and r. 



Appendix B. Inverse Compton 

In this appendix, we write down Boltzmann equations for inverse Compton process which 
determine the high energy photon spectrum. The electron energies before scattering and 
after scattering are given by E^ and E'^. E^ is used for the energy of scattered photon 
and for background photon. 

The diffuse extragalactic 7 ray flux has been observed by the COMPTEL and 
EGRET measurements. COMPTEL and EGRET observed 7 ray energy ranges from 
O.SMeV to lOOGeV. Electrons and positrons which scatter up background photons to 
this energy range should be highly relativistic. Thus, the number of collisions per unit 
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time per photon energy through inverse Compton process is written as 



X 



1 (Tq? \ 
2glng + (1 + 2g)(l - q) + 2^:^(1 " l)) ' (^-l) 



where is classical electron radius, F = 4e^E(./ml and q = E^/T{E(. — E^). The 
maximum photon energy is given by Ef.r/{1 + F). 

The Boltzmann equations for inverse Compton process is given by 



dt 



I, , dEJ,{E,) de,Me,)--r{e,,E,,E,) (B.2) 



f^(E'^) l-Ema^ f°° (fN , 

^gf^ = J dEefe{Ee) J ^^^^''^^^^ -^7 = -^e + £7 " E'^, -E'e) 

-fe{K) dt,h{t,) — {t,,K) (B.3) 



where -Emax = TTie^{E'^ — e^)/(4e^^ + rrie^ — 4e^-Ee) when the denominator is positive 
and Emax = oo when negative. The distribution function of the background photons at 
temperature is represented by /^(e^), 

vr^ exp(e^/r^) — 1 

To check the validity of these Boltzmann equations, we show the electron number 
conservation. 



r dE'^^I^ = r dEME,) [ ^ dE'^ r de,Me,)^-^{e„ E, + e,- E'^, E,) 

Jm^ ot Jm^ J E^j{i+r)+t-., Jo atah,^ 

- / dEME'J / de,Me,) — {e„E:) 

Jme Jo ^'^ 

= / dEefe{Ee) / dE^ / de^fb{e^) {e^,E^,Ee) 

Jme Je-y Jo atarj^ 



>=7 

j-co j-co .yy 

- j dE'ME'^) deMe,) — ie,,E'^) 

POO r-oo ^J^ 

= / dEefe{Ee) / de^fb{e^) — {e^,Ee) 

Jme Jo 

f°° f°° dN 

- / dE'ME'J / de,Me,) — {e„E'J 

Jme Jo 
= 

The energy conservation is also easily shown. 

dE.iE, - e,)^^n^ = j dEJ,iE,)j deMe,) 

rEeV/{l+V) 

X j dE^ {E^ - e^) ^^^^ (e^ , E^ , Ee) 
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Figure Bl. Photon spectra which are produced through inverse Compton process 
when high energy electrons are injected with = 10^ GeV and = 1 GeV^. From 
upper to lower the lines represent the photon spectra when 1 + z — I (thin solid line), 
10 (thin dotted line), 100 (thin dashed line) and 1000 (thick sohd line). Thick dotted 
line represents the initial electron spectrum. For simplicity, we have assumed here that 
background photon spectrum is monochromatic with — 2.7T^ 



In Fig. [Bl\ we plot the photon spectra which are produced through inverse Compton 
process when high energy electrons are injected. The photon spectra have quite a 
different form for different value of F. In the Thomson limit corresponding to F -C 1, 
the first two terms of the right hand side of Eq. (IB. ID is dominant. In this case, many 
low-energy photons are produced. In the extreme Klein-Nishina limit corresponding to 
F ^ 1, the last term of the right hand side of Eq. (IB. II) is dominant at larger q and 
photon spectra have peaks near the high-energy end |[53j. Our calculation is different 
from the steady-state method [531 El] • However, once electron spectrum is in a steady- 
state, our result is in good agreement with the result in [53l IM] . 
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